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in : 1 Introduction 

Using the full Belle T(5S) data sample of 121 fb _1 we have measured exclusive branch- 
ing fractions for the decays — > dW+Za*) - [lj. Assuming these decay modes sat- 
urate decays to CP-even final states, the branching fraction determines the relative 



Q ■ width difference between the CP-odd and -even eigenstates of the B s . 

<N ! 2 Event Selection 

> 

00 ' At the T(55*) resonance, the total number of B S B pairs produced is N = (7.11 ± 

1.30) x 10 6 . Three production modes are kinematically allowed: B S B S , B S B* (or 
B*B S ), and B*B*, where the B* decays via B* — > B s ^. In this analysis we obtain 
the signal yield for the last mode, for which the fraction (/ R »-5») is 0.870 ± 0.017 [2]. 

^ s & s 

We reconstruct B° s — > D* + D*~, D^DJ, and DfD~ decays in which D* + — >■ 
D+7 and P+ <j>(-> K + K-)ir+, K%(-> tt+tt") K + , K*°(^ K + tt~)K + , <f>(-> 
K + K-)p + (^ tt+tt ), Ifg(-). ii + tt-)K* + (^ K° s tt+), and K*°(-^ K + n~) i^* + (— )■ 

We select B s candidates using two quantities evaluated in the CM frame: the 



beam-energy-constrained mass M bc = y P^eam — Pb an d the energy difference AE = 
E B — E heam , where p B and E B are the reconstructed momentum and energy of the 
B s candidate and Pbeam is the beam energy. The fit region is selected to be M bc in 
[5.25,5.45] GeV/c 2 and AE in [-0.15, 0.10]_GeV. Because thej from B* -> P s7 is 
not reconstructed, the modes T(5S f ) — > B S B S , B S B* and B*B* are well-separated 
in M bc and AE. We expect only small contributions from B S B S and B S B* events 
and fix these contributions relative to B*B* according to our measurement using 
B° s D~7t + decays 0. 

Approximately half of selected events contain multiple candidates, typically due 
to random photons incorrectly assigned as D* s daughter photons. For these events we 
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select the candidate that minimizes the quantity q+n) \j-^ D s — M Da )/cr 

Edj [AM - AM)/(x AM ] 2 | , where AM = M~+^ - M~+ and AM = Mp. - M^. 

The summations run over the two daughters and the iV (=0, 1, 2) D* + daughters 
of the B® candidate. 

Background from e + e~ = w, <i, s, c) continuum events are rejected using 

a Fisher discriminant based on a set of modified Fox- Wolfram moments [I]. The 
remaining background consists of T(5S) ->■ fiW^W -»■ D+X, T(55) RBX, and 
fi s -»■ D±j(2317)DW, 5 S -»■ D^(2460)DW, and 5 S DfDJir decays. The last 
three processes are estimated to be small using analogous B d — > D sJ D^ branching 
frations and considered when evaluating systematic uncertainty due to backgrounds. 

We measure signal yields by performing a two-dimensional extended unbinned 
maximum-likelihood fit to the M^-AE distributions. For each sample, we include 
probability density functions (PDFs) for signal and background. 

The signal PDFs have three components: correctly reconstructed (CR) decays; 
"wrong combination" (WC) decays in which a non-signal track or 7 is included in 
place of a true daughter track or 7; and "cross-feed" (CF) decays in which a D* T is 
not fully reconstructed or fake. In the former case, the 7 from D* + — > Df'j is lost and 
AE is shifted down by 100 — 150 MeV; this is called "CF-down." In the latter case, 
an extraneous 7 is included and AE is shifted up by a similar amount; this is called 
"CF-up." All signal shape parameters are taken from MC simulation and calibrated 
using £?° D^~tc + and B° ->■ D^ + D~ decays. The fractions of WC and CF-down 
events within the fit region are taken from MC. The fractions of CF-up events are 
floated as they are difficult to simulate accurately (i.e., many B® partial widths are 
unmeasured). As the CF-down fractions are fixed, the separate DfD~, D^DJ, and 
D* + D*~ samples must be fitted simultaneously. 

The projections of the fit are shown in Fig. [T]and the fitted signal yields and the 
branching fractions are listed in Table [TJ The systematic errors are listed in Table [2J 



Mode 


Y (events) 


e (xl0~ 4 ) 


B (%) 


S 


DtD- 


33.11™ 


4.72 


0.581°,;^ ±0.13 


11.5 


D*±DJ 


44.511;! 


2.08 


1.8 ±0.2 ±0.4 


10.1 


d* s d* s 


24.4l|i 


1.01 


2.0 ±0.3 ±0.5 


7.8 


Sum 


102.01H 




4.3 ±0.4 ±1.0 





Table 1: Signal yield (Y), efficiency including intermediate branching fractions (e), 
branching fraction (B), and signal significance (S) including systematic uncertainty. 
The first errors listed are statistical and the second are systematic. 
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D* S D S 






Source 


+ CT 


—o 


+a 


—a 


+(T 


—o 


Signal PDF shape 


2.7 


2.2 


2.2 


2.4 


5.1 


3.8 


Bckgrnd PDF shape 


1.5 


1.3 


1.3 


1.4 


2.9 


2.8 


WC + CF fraction 


0.5 


0.5 


4.7 


4.5 


11.0 


9.7 


1Z requirement 


3.1 


0.0 


0.0 


2.7 


0.0 


2.1 


Best cand. selection 


5.5 


0.0 


1.5 


0.0 


1.5 


0.0 


7r ± /ir ± identif. 


7.0 


7.0 


7.0 


7.0 


7.0 


7.0 


Ks reconstruction 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


7T° reconstruction 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


7 






3.8 


3.8 


7.6 


7.6 


tracking 


2.2 


2.2 


2.2 


2.2 


2.2 


2.2 


Polarization 


0.0 


0.0 


0.8 


2.8 


0.6 


0.2 


MC statistics for e 


0.2 


0.2 


0.4 


0.4 


0.5 


0.5 


B{vf) 


8.6 


8.6 


8.6 


8.6 


8.7 


8.7 


N 

f -* 






18.3 
2.0 






Total 


22.7 


21.8 


22.7 


22.9 


26.2 


25.5 



Table 2: Systematic errors (%). Those listed in the top section affect the signal yield 
and thus the signal significance. 



3 Ar s Estimation 

In the heavy quark limit with (m b — 2m c ) — > and N c — > oo, the dominant contribu- 
tion to the decay width comes from B® — > £)(*)-£)(*)+ decays [5J |6]. Assuming negli- 
gible CP violation, the branching fraction is related to Ar s as Ar s /r s = 2B/(1 — B). 
Inserting the total B from Table [T] gives 

— - = 0.090 ±0.009 ±0.023, (1) 

s 

where the first error is statistical and the second is systematic. The precision of 
this result is similar to that of other recent measurements [7J [8] and consistent with 
theoretical predictions [10J. Two main uncertainties are the unknown CP-odd com- 
ponent in B° — > D* + D*~ decay and the size of contributions from three-body final 
states. The former is estimated to be only 6%. However Ref. [9] calculates significant 
contribution from B° s — > D^*' £)W J{{*) decays. This calculation predicts Ar s /T s from 
alone to be 0.102 ± 0.030, which agrees well with our result. 
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Figure 1: and AE projections of the fit result. The columns correspond to 

B° s D-D+ (left), B° s D^D^f (middle), and £° D*"!)^ (right). The red 
dashed curves show CR+WC signal, the blue and purple solid curves show CF, the 
grey solid curves show background, and the black solid curves show the total. 



4 Polarization Measurement of D* D* + 

We have also made the first measurement of the longitudinal polarization fraction 
(f L ) of D*~D* + . We select events using the same criteria as before, however we use 
a narrower range of M& c and AE (2.5cr in resolution) in order to minimize — » 
D* ± D* =F cross-feed. For these events we perform an unbinned maximum-likelihood 
fit to the helicity angles 9 l and 9 2 , which are the angles between the daughter 7 
momentum and the opposite of the B° s momentum in the D*^ rest frame. The 
angular distribution is [\A + \ 2 + |A_| 2 ) (cos 2 6\ + 1) (cos 2 6» 2 + l) + |A | 2 4sin 2 6' 1 sin 2 6 2 , 
where \A + \, \A_\, and \A \ are the helicity amplitudes. The fraction f L is given by 
IA)IV(IA)| 2 + l^+l 2 + l^-H- We measure 

h = 0.06 ±8i? ±0.03, (2) 

where the first error is statistical and the second is systematic. The systematic errors 
are from signal PDF shapes (+0.008, -0.010), background PDF shape (+0.007, -0.004), 
fixed WC fractions (+0.013,-0.015), fixed background level (±0.022), continuum 
suppression (+0.011), possible fit bias (—0.011), and MC efficiency due to statistics 
(±0.0004). The helicity angle distributions and fit projections are shown in Fig. [2j 



4 



-1 -0.5 0.5 1 -1 -0.5 0.5 1 

cos0! cos9 2 

Figure 2: Helicity angle distributions and projections of the fit result. The green 
dashed (red dashed-dotted) curves show the transverse (longitudinal) components, 
the blue dotted curve shows background, and the black solid curve shows the total. 



In summary, using the branching fractions of B® — >■ £)(*)+ decays, we mea- 
sured Ar s /r s assuming no CP violation, among other theoretical assumptions. We 
have also made the first measurement of the B° s — > D*~D* S + longitudinal polarization 
fraction. 
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